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Abstract 

The Lorenz attractor is one of the best known examples of applied mathematics. However, 
much of what is known about it is a result of numerical calculations and not of mathematical 
analysis. As a step toward mathematical analysis, we allow the time variable in the three 
dimensional Lorenz system to be complex, hoping that solutions that have resisted analysis on 
the real line will give up their secrets in the complex plane. Knowledge of singularities being 
fundamental to any investigation in the complex plane, we build upon earlier work and give a 
complete and consistent formal development of complex singularities of the Lorenz system using 
psi series. The psi series contain two undetermined constants. In addition, the location of the 
singularity is undetermined as a consequence of the autonomous nature of the Lorenz system. 
We prove that the psi series converge, using a technique that is simpler and more powerful than 
that of Hillc, thus implying a two-parameter family of singular solutions of the Lorenz system. 
We pose three questions, answers to which may bring us closer to understanding the connection 
of complex singularities to Lorenz dynamics. 
Keywords: Lorenz attractor, psi series, complex singularities. 
AMS: 34M35, 37D45. 

1 Introduction 

The nonlinear system of equations 

d y oq 

— = Zax — y — xz 
dt 

dz 

— = -8 2 /3 + xy, (1.1) 

which is named after Lorenz, gives the best known example of a strange attractor. Lorenz |21[ [22] 
derived this system to argue that the unpredictability of weather is due to the nature of the solutions 
of the Navier-Stokes equations and not due to stochastic terms of unknown origin, his point being 
that a deterministic system could possess an attracting and invariant set on which the dynamics 
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is bounded and linearly unstable. When such strange attractors exist, trajectories are chaotic and 
appear random. 

While Lorenz |21j p. 141, 1963] could write that the atmosphere was not normally regarded as 
deterministic, we now know that the incompressible Navier-Stokes equations by themselves explain 
a remarkable wealth of turbulence phenomena including coherent motions in the near-wall region, 
the law of the wall, intermittency, and vortex structures in fully developed turbulence [2]. The 
density and temperature of the atmosphere vary with altitude, and there is significant electrical 
activity in the atmosphere that is sustained by about 40,000 thunderstorms that occur around the 
world in any single day [61 Chapter 9] . If we nevertheless think that the physics of the atmosphere 
is deterministic, Lorenz and his system are partly responsible. 

Lorenz's point of view was dynamical. He viewed the state of (jl.ip as a point in i? 3 and its 
solutions as trajectories in R 5 . The dynamical point of view has overwhelmingly dominated work 
on the Lorenz system and Lorenz's original paper [21] has remained an outstanding introduction to 
dynamics. In it, a careful reader can find discussions of numerical errors, of concepts of stability, of 
symbolic dynamics (aspects of which Lorenz seems to have rediscovered for himself), of the density 
of periodic solutions on the Lorenz attractor, and of the fractal nature of the Lorenz attractor. 

The point of view in this paper, unlike Lorenz's, will be mainly function theoretic. We view t in 
(jl.ip as a complex variable and x, y, z as analytic functions of a complex variable. Our interest is 
in triples of analytic functions which satisfy (II. ip . Our hope is that an investigation in the complex 
plane will open a route to the mathematical analysis of the Lorenz system. 

For the most part, we deal with certain singular solutions of the Lorenz system, which will be 
introduced momentarily. As the right hand side of the Lorenz system (jl.ip is analytic, every solution 
of the Lorenz system admits analytic continuation to the complex plane. For some solutions, the 
analytic continuations have singularities of the form we deal with, as indicated by numerical results 
summarized in Section 5. In the second part of this introduction, we pose three questions to help 
connect the complex singularities with Lorenz dynamics. 

From residue integration, the method of steepest descent, and the use of deformation of contours 
to effect analytic continuation of certain special functions, we know that knowledge of singularities 
is often useful to investigations in the complex plane. This observation explains our focus on 
singular solutions of the Lorenz system. 

1.1 Psi series solutions of the Lorenz system 

The most common types of singularities are poles, algebraic branch points, and logarithmic branch 
points. The singularities of the Lorenz system that we examine are of none these types, but are 
given by psi series representations. 

Definition 1.1. A logarithmic psi series centered at to is a series of the form X^n^-TV-P"! 7 ?)^ - to) n > 
where N is an integer, n = log(b(i — to)) an d each p n is a polynomial in n. In the definition of n, 
b is a complex number with \b\ = 1, with b = ±i often being convenient choices. 

Throughout this paper, log will denote the principal branch of log. The choice of the branch 
is ultimately immaterial but taking n = log(— i(t — to)) instead of r\ = log(i — to) leads to more 
convenient branch cuts if 9(to) < 0, as we explain in Section 3. For a slightly different definition of 
logarithmic psi series, along with definitions of psi series of other types, see [131 Chapter 7.1]. The 
only type of psi series that arises in this paper is the type given by Definition II .lj and by psi series 
we refer to that definition only. 
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The psi series of Definition 11.11 are like the Laurent series, except that the coefficients are 
polynomials in 77 instead of being constants. For that reason, the psi series singularities were called 
pseudopoles by Hille Even though the coefficients are polynomials in 77, each nonzero term of 
the logarithmic psi series dominates the following term in magnitude in the limit t — ► t$. 

In an intriguing and original pair of papers, Tabor and Weiss [32] and Levine and Tabor [19] 
considered psi series solutions of the Lorenz system (II. ip . The psi series they used were expressed 
as a double sum. Below we give the psi series in a different form: 

X(t) = + + P lW(* " *0) + P ^){t ~ tof + ■■■ 

t — 1 

v® = + Qofa) + Qi(v)(t - to) + Q2(v)(t - 1 ) 2 + ■■■ 

z ^ = jri$ + tt^ + Ro ^ + ~ to) + R ^ ~ to)2 + • • ■ (L2) 

Here the Pi, Qi, and Ri are polynomials in 77 where rj = log(6(t — to)) as i n Definition 11.11 As 
the Lorenz system is autonomous, to is an arbitrary complex number. The fact that the leading 
powers of (t — to) in the three series in (jl.21) are —1, —2, and —2 may be guessed by substituting 
poles (t — to) - ") (t — to)~P, (t — to) -7 for x, y, z into the Lorenz system and then solving for a, (3, 
7 by matching the order of the left and right hand sides [32j . This test-power method |13|. p. 90] 
does not always work and can be tricked into failing for the Lorenz system with a linear change of 
variables. 

Melkonian and Zypchen [21] have recast the psi series of Tabor and Weiss [32] into the formalism 
of Hille [UJ. The formal development of psi series that we give in Section 3 is similar to that of 
Melkonian and Zypchen [21], but improves that of Melkonian and Zypchen in two respects. Firstly, 
the development in Section 3 shows the dependence on undetermined constants C and D explicitly, 
pointing out the occurrence of rj and C in the group (77 + C). Secondly, we prove that the degrees of 
P-m+i, Qm, Rm are given by L^^^J for m = 0, 1, . . .. The proof hinges on a surprising cancellation 
for m = 2. It is important to get such details fully right if a mathematical theory is to be set up. 
As Hille |13, p. 68] pointed out, "constants of integration play a remarkable role in the advanced 
theory of nonlinear DEs." In addition, a complete formal calculation is essential for a fully correct 
convergence proof. 

The first few coefficients of the psi series (|1.2p are listed in Table [TJ It is evident that 77 and C 
always occur in the group (77 + C). If D were real, the coefficients of the polynomials in (77 + C) 
listed in that table would all be either pure imaginary or real. 

The following is one of our main theorems. It reappears in a more specific form in Section 4, 
where it is proved. 

Theorem 1.1. The psi series (11, 2D . some of whose coefficients are listed in Tabled satisfy the 
Lorenz system (jl.ip in the disc \t — to| < r for some r > 0, but with the singular point t = to and 
a branch cut deleted from the disc. The constants C and D are undetermined. 

The proof of this theorem is valid for any choice of the undetermined constants C and D, but 
the estimate for r depends upon the choice. A key step in its proof is to show the convergence of 
the psi series. 

An important aspect of the convergence of the Lorenz psi series is not brought out in Theorem 
11.11 As evident from the appearance of 77 in Definition ll.il a typical psi series will have logarithmic 
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Table 1: Coefficients of the psi series of (jl.2p . with 77 as in Definition 11,11 Evidently, the degrees 
of Pm+i, Qm, Rm are L^^r^J f° r m = 0,1,2,3. The other valid choice for coefficients of (11. 2j) 
is obtained by changing the signs of all the PjS and QiS, while leaving the RiS unchanged. The 
constants C and D are both undetermined. 

branch points in the t-plane. To get around the multiple- valuedness, Theorem 11.11 fixes a branch 
cut in the i-plane. The branch cut can be dispensed with by parametrizing the Riemann surface 
using i]. A discussion of convergence in the ry-plane is found in Section 4 (see Figure[3]in particular). 

Hille's "frontal attack" to prove convergence of psi series can be modified to apply to the Lorenz 
system 124] . In an appendix, Hille [12] pointed out that his technique could only handle the 
Emden-Fowler system (see Section 2) with p = 2, while a more complicated technique due to Smith 

[30] could handle p = 2,3, The technique we use in Section 4.1 is also a frontal attack, but it is 

a good deal more transparent than Hille's approach. In place of an elaborate analytic set up and 
an inductive hypothesis to bound the coefficients of the psi series, we use the Laplace transform, 
elementary combinatorics, and an elementary implicit function theorem. Our technique seems to 
extend to all the cases handled by Smith [30J. Detailed comments on this point are found in Section 
4.2. 

1.2 Complex singularities and Lorenz dynamics: three questions 

From Theorem 11.11 we get a two-parameter family of singular solutions of the Lorenz system (jl.ip . 
The form of the singular solutions is given by the psi series (jl.2p and the two undetermined constants 
C and D are shown in Table [TJ The location t$ of the singularity can be anywhere in the complex 
t-plane. 

For some definite integrals, the singularities of the integrand and Cauchy's residue theorem 
imply the value of the integral. So we ask, what do the singular solutions of the Lorenz system 
tell us about the dynamics in R? for real time? As the analytic theory of solutions of the Lorenz 
system is still in its infancy, a complete answer to the question cannot be given. Nevertheless, the 
question merits a thorough discussion. 

Many beautiful visualizations of the Lorenz attractor are found on the INTERNET. The visu- 
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Figure 1: The periodic orbit in the first plot is labeled AB to indicate the sequence in which it 
moves between the A quadrant (with x < —16.432, y < —16.432, z = 27) and the B quadrant (with 
x > 16.432, y > 16.432, z = 27). Each filled circle is directly below a singularity in the complex 
t-plane. In the middle are plots of x(t) (solid), y{t) (dashed), z{t) (dotted) against real t. In the 
rightmost plot, the location of the complex singularities of AB that are closest to the real line are 
marked as crosses. The orbit AB is computed with 547 digits of precision. 

alizations originally offered by Lorenz [21] are packed with information and are models of concision. 
The Lorenz attractor is a butterfly-like subset of R 3 . Except for the fixed points, all trajectories 
either approach the attractor as t — ► oo or are already on it. 

Figured] shows the periodic orbit labeled AB, which resides on the attractor. A great advantage 
of computing such orbits, as opposed to arbitrary trajectories, is that the computations take on 
a definite character that makes it possible to report them precisely. As already mentioned at the 
beginning of this introduction, periodic orbits are believed to be dense in the Lorenz attractor. Such 
orbits can be computed with great precision. The locations of the complex singularities shown in 
the rightmost plot of Figure [U were obtained by computing the orbit AB with more than 500 digits 
of precision. 

A worthy goal for the analytic theory of the Lorenz system is a proof of existence of periodic 
solutions (x(t),y(t), z(t)) of the Lorenz system (jl.ip . where we seek a proof that is based solely on 
mental conceptions. There is a definiteness to seeking periodic solutions as already pointed out. 
In addition, periodic orbits are key to extracting order from chaos, to borrow an expression from 
Strogatz [21]. For instance, Figure which illustrates the fractal property of the Lorenz attractor, 
was obtained by computing periodic orbits. The plots were computed in parallel on a machine with 
two quadcore 2.33 GHz Xeon processors. The plots took a day or two of computing. For the theory 
behind such computations, see [33] and [34J. 

A proof of existence of periodic solutions of the Lorenz system (jl.ip appears to be far away. 
We formulate three questions to serve as more immediate goals for the development of the analytic 
theory of the Lorenz system. 

Question 1.1. Are all singular solutions of the Lorenz system given by psi series expansions (II. 2j) 
with suitable choice of the undetermined constants C and D ? 

The role of the undetermined constants C and D is partly shown in Table [TJ Their role in the psi 
series is clarified further in Sections 3 and 4. Lorenz [21J gave arguments that partially imply that 
a real solution of the Lorenz system cannot become singular in finite time. The implication covers 
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Figure 2: Fractal property of the Lorenz attractor. (a): The intersection of an arbitrary trajectory 
on the Lorenz attractor with the section z = 27. The plot shows a rectangle in the x-y plane. All 
later plots ((b) and above) zoom in on a tiny region (too small to be seen by the unaided eye) at 
the center of the red rectangle of the preceding plot to show that what appears to be a line is in 
fact not a line. These plots, and the plots of |33|. l34j of which these plots are a refinement, appear 
to be the only plots made of the fractal structure of the Lorenz attractor. 



both increasing and decreasing time. In Section 5, we give a complete proof of that result. Thus for 
solutions of the Lorenz system that are real for real t, the locations to of the complex singularities 
must have a nonzero imaginary part. In fact, Foias and others [HJ Theorem 2.3] have proved that 
for solutions on the Lorenz attractor the imaginary part of the location of the singularity in the 
complex t-plane must exceed 0.037 in magnitude. For an investigation of the backward in time 
behavior of the Lorenz system (for real data), see the paper by Foias and Jolly [7]. 

The techniques used to deduce psi series solutions of the Lorenz system are not of much use for 
answering Question 1 1.1 1 However, if to is an D singular point of the Lorenz system, then \x(t)\ + 
\y(t)\ + \z(t)\ — > oo as t — > to, as implied by a slightly stronger theorem proved in Section 5. 

For analytic functions such as the gamma and zeta functions, analytic continuation into the 
complex plane is an important step in understanding the true nature of those functions [25]. The 
question of analytic continuation is important in the theory of differential equations in the complex 
plane as well |13j . These observations motivate us to ask the following question. 

Question 1.2. Do solutions of the initial value problem for the Lorenz system with (x(0),y(0), z(0)) 
being finite (but possibly complex) admit of analytic continuation to the entire complex t-plane except 
for branch points? 

An affirmative answer to Question 11.11 appears to imply an affirmative answer to Question ll.2i 
The process of analytic continuation can be blocked by singularities. But if all singularities are 
given by psi series of the form (|1,2[) , Theorem 11.11 implies that we can continue around any such 
singularity into a disc of finite radius around that singularity (radius is r in the theorem). The 
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possibility where a succession of psi series singular solutions of decreasing radii of convergence 
accumulate on another singular point is easily ruled out, if the answer to Question 11.11 is yes. 

Singular solutions given by psi series representations exist for plane quadratic systems as well 
as plane polynomial systems [121 130] . Such planar systems certainly cannot exhibit chaos [31J. 
The dynamics of planar systems is tightly circumscribed by results such as the Poincare-Bendixson 
theorem. Unlike the Lorenz system, the planar systems considered by Hille |12j and Smith [30] can 
have real solutions that develop singularities in finite time. Yet one is probably justified in thinking 
the mere existence of singular solutions represented by psi series is unlikely to tell us anything 
about the chaotic nature of the Lorenz system. 

This is perhaps the place to comment on the three free parameters with which the Lorenz 
system is usually written, but which are given the values used by Lorenz [21] in (II. ip . The three 
parameters correspond to the Rayleigh number, the Prandtl number, and the system size for the 
convection PDE from which the Lorenz system was derived. With regard to the choice of these 
parameters, there are three cases for which the Lorenz system admits a Laurent series as a solution 
[29l [32] . There are five other cases, due to Segur [29] and Kus [15] . for which time-dependent 
integrals of motion are known. In their pioneering work, Tabor and Weiss [32] considered the 
connection between integrability and the type of the singularities. For another discussion of the 
connection between psi series and integrability, see [3]. 

In addition to the integrable cases, there are a number of other regions in parameter space 
where the Lorenz system has non-chaotic dynamics yet admits singular solutions with psi series 
representation. In these instances, it is quite possible that even though the real-valued dynamics is 
non-chaotic, more varied solutions exist when complex numbers are allowed. In the case of plane 
polynomial systems, although the differential equations cannot have chaotic solutions that are real 
[T| 131]. the equations may have chaotic solutions that are complex. 

It is not entirely clear how the nature of the singularity can be connected to chaotic dynamics. 
It is perhaps significant that only real solutions have a bearing on dynamics. Therefore we ask the 
following question. 

Question 1.3. If a psi series solution of the Lorenz system (jl.ip of the form ([ 1 . 2 j) is obtained by 
analytic continuation of a solution that is real for real t, what constraints must C, D and to satisfy? 

The detailed development of psi series found in Section 3 and partly shown in Table Q] could 
help answer this question. Numerical computations are also likely to be useful. A suspicion of ours 
is that the undetermined constant D is real for the psi series singularities of Question 11.31 

2 A brief history of early work on psi series 

The equation of Briot and Bouquet 



where F is a polynomial with quadratic and higher terms, seems to be the simplest differential 
equation whose singularities are given by psi series. Dulac [4, p. 368, 1912] and Malmquist [231 P- 
19, 1921] (also see Theorem 11.3.1 of [13]) proved that the general solution of (12, ip around t = 
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is given by a convergent psi series if p is a positive integer. For generalizations to higher order 
Briot-Bouquet equations, see [17] . 

In the last decade of his life, Einar Hille |101 111] \12\ 113] interested himself in the Emden- Fowler 
equation cPy/dt 2 = t~ 2 / p y l+2 / p with p > 1 being a positive integer. The Emden- Fowler equation 
originally arose in cosmology. The special case p = 2 is the Thomas-Fermi equation, which arose in 
atomic physics. After sixty years of encounters with differential equations, Hille wrote a splendid 
book on ordinary differential equations in the complex plane \13\ 1976]. The last chapter of that 
book gives an outline of the work of Hille and Russell A. Smith [30J on psi series singularities of 
the Emden-Fowler equation. The techniques involved are highly relevant to the Lorenz system. In 
Section 4, we point out that some of the theorems of Hille and Smith admit simpler proofs using 
an approach introduced in that section. 

From Hille's illuminating bibliographic discussions [13], it is clear that Dulac [5j 1934] was a 
central figure with regard to psi series, with Horn [14} 1905] being another early contributor. Hille 
does not mention Dulac's claim about one of the Hilbert problems, however, and indeed that claim 
was mistaken [15] . It appears that the error was related to a subtlety in the interpretation of psi 
series in the complex plane [15] . 

3 Formal development 

The formal development of psi series has a history that goes back a hundred years or more. All 
formal developments proceed in a similar way — one begins with psi series and then determines their 
coefficients using a recursion. In two of his papers, Hille [111 [T2] gave clear and detailed formal 
developments. Our derivation is quite similar, but is more careful about subtleties such as the 
choice of the branch of log, the degrees of the polynomials Pi, Qi and Ri in (11 .2[> . and the role of 
the undetermined constants (C and D in Table [1]). 

Since the Lorenz system (jl.ip is autonomous, the choice of the location to of the singularity is 
arbitrary. For the sake of definiteness and because the primary interest is in solutions that are real 
for real t, we may assume 9ft(io) < 0. and take rj = log(— i(t — to)) to obtain a branch cut that does 
not intersect the real axis. However, nothing changes if to is arbitrary and some other branch cut 
is chosen for defining r/. The choice of branch cut is equivalent to the choice of b in Definition 11.11 

The form of the singularity is assumed to be given by (jl.2p : 

oo oo oo 

X(t)= Y, Pm(v)(t-t ) m V(t)= Yl Qm(v)(t-t ) m Z(t)= Y, Rm(v)(t-t ) m , (3.1) 
m=— 1 m=— 2 m=— 2 

where the P m , Q m and R m are polynomials in 77. We arrived at this form based on numerical work 
summarized in Section 5. However, the credit for discovering the form of the psi series singularities 
of the Lorenz system belongs for the most part to Tabor and Weiss [32] , 
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Substituting (|3.ip into (|l.ip and denoting derivatives with respect to rj by a prime, we get 

oo oo 

Y, (KM + mP m (rj))(t - hr- 1 = 10Q_a(* - ^o)" 2 + E ( W Qm(v) - ±0P m ( V ))(t - t ) m (3.2a) 

rn— — 1 m=— 1 

oo oo oo 

£ (QU^ + mg^x*- to)" 1 " 1 = 28 J2 Pm(v)(t-t ) m - E Qm{n){t - h) m 



m=-2 

oo / m+2 



E E Pj{l)R<n-M (*-M m (3-26) 

m=-3 \j=— 1 / 



8 x 



E (^(r ? )+m J R m (ry))(t-t ) m - 1 = -3 E Rm(v)(t-t ) m 

m— — 2 m——2 

oo / m+2 \ 

+ E E Pi(v)Qm-Av) (t - to) m (3.2c) 

m=-3 \j=-l / 

For the psi series on either side of (3.2), a nonzero term with m = mi is greater in magnitude 
than an m = mz term in the limit t — * to if mi < m<i- Therefore it is formally consistent to equate 
powers of (t — to) in increasing order. 

Equating coefficients of (t — io)~ 2 in (|3.2aj) and of (t — to) -3 in (|3.26j) and (|3.2cp . we get 
P'_ x - P-i = WQ-2, Q'- 2 - 2 Q-2 = -P-1R-2, and R'_ 2 - 2P- 2 = P-xQ-2- The degree of P_i 
and Q-2 in r\ must be the same, while the degree of R-2 must be twice that degree and the degree 
of Q-2 must be the sum of the degrees of the other two. The only possibility is for all the degrees 
to be zero. We get 

(P_i,Q-2 ) fl-2) = (2» ) -*/5,-l/5) or (-2^/5,-1/5). (3.3) 

We consider only the first possibility for now, but will account for the second possibility in Lemma 
1531 

The next set of equations is Pq = 10(Q_i — P-i), Q'-\ = Q-i — %R-\ + Po/5 + i/5, and 
R'_i = R-i + 2iQ_i — iPo/5 + 8/15. The only solution polynomial in r/ is given by 

(P , Q-i, R-i) = (71t/9, 2i, 17/9). (3.4) 

For m = 0, 1, 2, . . ., we equate powers of (t — to) m in (|3.2a|) and powers of (t — to) m ~ l in ()3.2 &[) 

X' m = A m X m + F m (r]), (3.5) 



The eigenvalues of A m are —m + 2, — m, and — m — 3. If the linear system (|3.5p is diagonalized using 
the eigenvectors of A m as a basis, it turns into three scalar equations of the form d^/dr] = a£+f(r]), 
with q being — m + 2 or — m or — m — 3 and with / being a polynomial in each case. If a ^ 0, we 
have a unique polynomial solution for £(77) whose degree is the same as that of /. 



and (|3.2c|) 


to get, 




where 








'Pn+l\ 






Qm 1 , 


j4rn — 1 
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We can have a = if and only if m = or m = 2. Thus if m / and m / 2, we can assert 
that (13. 5p has a unique polynomial solution X m and the degree of that solution in rj is the same as 
that of F m . 

In the case m = 0, -F m is a constant and the three scalar equations are of the form d£/dr] = 
2£ + /3i, d^/dr] = — 3£ + /?2 and d^/dr/ = /3$, where the are known constants. The only admissible 
solution of either of the first two equations is a constant. The last equation however has the solution 
+ C), where C is an undetermined constant. If the eigenvectors of Aq are multiplied by the 
respective solutions and summed, we get 

-9880i/8l\ / \ 

-988i/81 (rj + C)+ -349i/81 , (3.7) 
/ V 1385/54 / 




where the factor multiplying (rj + C) is the eigenvector of Aq that corresponds to the eigenvalue 
— m = 0. 

The matrix A m has a zero eigenvalue again when m = 2. In this case, the degree of F m in rj is 
2. We would expect the polynomial solution X m of (|3.6jl to be cubic. However, the component of 
F m along the eigenvector of A m corresponding to the eigenvalue — m + 2 = is zero (with regard 
to this point compare (2.9) of [32]). Therefore P3, Q2 and R2, which make up X2, are all quadratic 
in rj as shown in Table [TJ A new undetermined constant D enters at this stage. If P3, Q2 and R2 
were cubic and not quadratic, L 21 ^] i n the lemma below would be replaced by L 21 ^] + L^^J • 

Lemma 3.1. The degrees of the polynomials P m +i(r]), Qm(v) an d Rmiv) are a ^ m °st [ m ^ 2 j for 
m = 0,1,2, .... 

Proof. For m = 0, 1, 2, the lemma can be verified explicitly using Table [TJ If the maximum degree 
of a component of is dk for < k < m, (|3.4p and (|3.6p imply that the degree of F m is at most 

max (dj-i + dm-j-l), 

0<j<m 

where we assume m > 3 and take d_\ = 0. We use the inductive hypothesis and note 

dj-i + dm-j-i < + L 2 J - ^~^T~ * 

for < j < m to complete the proof. The second inequality above is an equality for odd j. □ 

It appears as if the degrees in Lemma 13. II are actually equal to L 222 ^]- To prove as much, one 
has to rule out cancellations that can happen in a variety of ways, which may or may not be worth 
the trouble. Below we give a formula for the polynomial solution X m of (|3.5p that is easily derived 
using the variation of constants formula and integration by parts: 

1 m + 2 1 

2 M TP 

j=0 1 

for m > 3. The correctness of (|3.8p can be verified by direct substitution into (|3.5p . 
The lemma below summarizes the discussion in this section. 
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Figure 3: Schematic plot of the location of the singularities in the t-plane for an orbit such as AB. 
The singularities are shown as red spots and the branch cuts are dashed. Only singularities within 
a single period are shown in the t-plane (compare Figure H]). 



Lemma 3.2. (i) For the coefficients P m +i, Qm> R-m shown in Tabled for —2 < m < 3 and 
defined for m > 3 by (|3.6p and f|3.8j) . the psi series (|3.ip (or (|1.2p ) satisfy the Lorenz system 
(jl.ip formally. The location of the singularity t$ is arbitrary and two undetermined constants, 
C and D, occur in the psi series. The constant C and n always occur in the group (n + C). 

(ii) Another formal solution is obtained by flipping the signs of all the Ps and the Qs while leaving 
the Rs unchanged. 

(Hi) For the solution to be formally valid, n can be defined as log(6(t — to)) for any complex number 
b with \b\ = 1. 

Proof. For the part about flipping signs, note that the Lorenz system is unchanged by the trans- 
formation (x,y,z) — » (— x, — y, z). More specifically note that flipping signs of the Ps and the Qs 
changes the sign of the first two components of F m in (13, 6p but not that of the third component. 
This other formal solution accounts for the second possibility in (|3.3[) . 



If the psi series singularity is an analytic continuation of a solution that is real for real t, the 
location to of the singularity must be off the real line (see Section 5). According as 9(to) < or 
Sj(to) > 0, the choices b = — i or b = i give branch cuts that do not intersect the real line, as shown 
in Figure El 

4 Proof of convergence 

Hille's [TT] proof of the convergence of psi series solutions relies on the formula 



for the solution X m of (|3.5p which is polynomial in n. A similar formula is fundamental to the 
approximation of strange attractors, including Lorenz's, by algebraic sets in the work of Foias, 
Temam and others [8j E] . 



□ 
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Our proof of convergence does not use Hille's formula, but instead relies on the Laplace transform 
and other devices. In the second part of this section, we remark that our technique will likely give 
simpler proofs for certain theorems of Hille and Smith. In one instance, our technique can probably 
be used to prove a theorem that has been stated but not proved completely. 

4.1 Psi series solutions of the Lorenz system 

If p is a polynomial in (77 + C), we define \p\ as the sum of the absolute values of its coefficients. 
Since 77 and C always occur in the group (77 + C), we can think of C as being subsumed by 77. 
\X m \ is defined as the maximum of |P m+ i|, \Q m \ and \R m \. For m > 2, \X m \ will depend upon 
the undetermined constant D. The key to the proof of convergence of the psi series (13. 1|) is a 
bound of the form \X m \ < KxK™, where K\ and K2 are positive constants that depend upon the 
undetermined parameter D. 

For F m defined by (|3.6p . \F m \ is the maximum of |-| over its three components, each of which is 
a polynomial in (77 + C). We begin with the following easy lemma. 

Lemma 4.1. For m > 3, 

m— 1 

\F m \ < 30 |X m _i| + 28 |X m _2| + |JT m _j_i| l-Xj-il . 



Proof. If p and q are polynomials in 77 + C, \pq\ < \p\ \q\ and \p + q\ < \p\ + \q\. Repeated use of 
those inequalities with the definition (|3.6p of X m and F m gives 

m 

\Fm\ < 10 |X m _i| + 28 |X m _2| + \X m -j-i\ |-Xj-i| • 

i=o 

The lemma results when the j = and j = m terms are moved out of the summation while using 
Tabled] to note that |X_ X | < 10. □ 

For matters related to the existence and uniqueness of the Laplace transform that arise implicitly 
in the proof below, see |35j . In the lemma below, we only treat polynomials in 77 (assuming C = 0), 
but the lemma still applies when 77 and C occur in the group (77 + C) and 0. 

Lemma 4.2. Let a be a complex number with \a\ > 1 and let f(rj) be a polynomial in 77. Let £(77) 
be the polynomial solution of the differential equation 

^ = 0$ + /(!/). (4.1) 

If the polynomial f(r]) is of degree n, assume \a\ > a(n + 1/2) for some a > 1. Then 

KI^A-^I/I- ^ 

a a — 1 
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Proof. Let f(rj) = fo + fir] + ■ ■ ■ + f n V n - To take the Laplace transform of (|4.ip . we multiply (|4.ip 
by e _r?s and integrate from 7/ = to 77 = 00. We get 



s£(s) - a£(a) = £(0) + — + 



/o , , 2!/ 2 



+ 



+ ••• + 



n!/ w 
„n+l ' 



Rearranging, we have 



/o 



+ 



+ • • • + 



s — a ' (s — a)s ' (s — a)s 2 ' (s — a)s n+l 

All terms on the right hand side above except the first are rewritten using the identity 



1 



1 



(s — a)s k a k (s — a) a k s a k 1 s 2 



as 



k ■ 



In the resulting expression, £(0) is chosen to cancel all the l/(s — a) terms to get a polynomial 
solution. We then have 



1 1 

+ ^7 + '-- + 



1 



A:=0 
n+1 

■E 

k=l 



a k+l s a k s 2 



1 / (fc - !)!/*_! , klf, 



+ 



a 



a" 



as 



+ ■•• + 



fc+l 



n\f n 



a 



n+2-k 



(4.3) 
(4.4) 



The coefficients of £(77) are evident from inspecting the summations (|4.3p and (|4.4p . From the 
summation (|4.3h and the inverse Laplace transform, we get 



ki<£i/*i 



fc! 



+ 



A;! 



fc=0 



0!|q /c+1 | l!|a fc 



+ •■■ + 



fc! 



+ 



A;! 



(fc - 1)! |a 2 | fc! |a 



11 I „ I / k 
Z-^ \ a \ \ /—< j\ \ n k-j\ 



k\ 



k=0 



\j=0 



1. (4.5) 



To clarify the calculation that gives (|4.5p . let us consider the special case d^/dr] = a£ + rj k . Its 
unique polynomial solution is £ = rj k /a — krj k ~ 1 /a 2 — • • • — k\/a k+1 and this |£| corresponds to the 
fcth term in (|4.5p . 

Next we bound fc!/j! |a fc_ - 7 j for < fc < n and < j < fc. 



fc! 



j -A 



j-fc 



fc(fc — 1) . . . (j + 1) 

(fc(j + l))((fc - l)(j + 2))((fc - 2)(j + 3)) . . . L 



< 



k + j + 1 



k-j 



< 



fc + 1/2 



k-j 



< 



n+1/2 



k-j 



a 



a 



2 |ct| / 
< l/a k ~ j . 

In the second line above, the last factor L is either (fc + j + l)/2 or ((fc + j)(k + j + 2)/4). The first 
inequality in the third line is obtained by applying the inequality xy < ((x + y)/2) 2 repeatedly. 
The inequality in the last line uses the assumption \a\ > a(n + 1/2) made in the statement of the 
lemma. 
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Returning to (|4.5p . we have 



|£|<^(l + l/a + l/a 2 + ---), 




which completes the proof. □ 

The inequality in the lemma below is not strict mainly because \F m \ = is not ruled out. 
Lemma 4.3. For m>8, \X m \ < 192 \F m \ /(m - 2). 
Proof. We take the matrix of eigenvectors of A m defined in (|3.6h to be 

V - I -:]//2 / i 

where the columns are ordered to correspond to the eigenvalues — m + 2, —m, and —m — 3, respec- 
tively. 

If (13. 5p is rewritten using a similarity transformation that turns A m into a diagonal matrix, we 
get three scalar equations 

-r- = UiZi + /i, 

(277 

for i = 1,2,3, where (01,02,03) = (-m + 2, -m, -m - 3), (ft, f 2 , f 3 )' = ^ -1 ^m, and X m = 
V(£i, £2, £3)' (the prime denotes transpose). 

To apply Lemma B~2l to each of the scalar equations, we use Lemma lBTTI and take n = [(m+2)/2\ . 
In addition, we choose an a > 1 such that 

|oi| = \m - 2| > a(m + 3)/2 > a(n + 1/2). 

The choice a = 12/11 works for m > 8. Thus we get |&| < 12\fi\/\ai\ < 12|/j|/(m - 2) for 
i = 1,2,3. 

We have \fi\ < \\V~ 1 \\ 00 \F m \ for i = 1,2,3 and [X m | < HVHoo maxd^ , |£| 2 , |£| 3 ). Combining 
the inequalities, we get 

I X m J < — 1 1 V 1 1 oo 1 1 V 1 1 00 I Fm I • 

m — 2 

The proof is completed by verifying that ||y||oo||^~ Hoc = 16. □ 

The lemma below is crucial to showing that the psi series expansions which formally satisfy the 
Lorenz system by Lemma 1331 are convergent. Its proof is structured to be transparent, but does 
not give the best constants. 

Lemma 4.4. For positive constants K\ and K% which depend upon the undetermined constant D 
of LemmalEM \X m \ < KiK% 1 for m = 0, 1, 2, . . . 

Proof. By Lemmas 14.11 and 14. 31 we have 

30x192 28 x 192 192 ^ 

m — 2 m — 2 m — 2 ^-^ 

3=1 
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for m > 8. If we define x m = \X m \ for m = 0, 1, . . . , 7 and, for m > 8, 

m— 1 

x m = 960x m _i + 896x m _ 2 + 32 ^ (4.6) 

i=i 

then |X m | < x m (after noting 192/6 = 32 and so on). 

Let f(Z) = J2m=o x mZ m be the generating function of the x m sequence. Using (|4.6p . we get 

f(Z) - (c + c x Z + • • • + c 7 Z 7 ) = 960Z/(Z) + 896Z 2 /(Z) + 32Z 2 f{Z) 2 . (4.7) 

In (14. 7p . the constants Co . . . C7 account for the fact that the recurrence (|4.6j) is valid only for m > 8. 
They are put in to get xo, x%, . . . , X7 as the coefficients of Z°, Z 1 , . . . , Z 7 , respectively. They can 
be determined explicitly (compare Tabled]); for instance, cq = xq = \Xq\ and c\ = x\ — 960xo = 
|Xi| — 960 |JTo|. Because IX2I , . . . , IX7I depend upon D, so will C2, . . . ,cj. 

The implicit function theorem implies the existence of a unique analytic function with /(0) = xq 
that satisfies (14. 7h — if all terms of (14. 7p are moved to the left and 'f(Z)' is treated as a variable, 
the partial derivative of the left hand side with respect to / is 1 when Z = 0, thus verifying the 
derivative condition of the implicit function theorem. Therefore f(Z) is the generating function of 
the x m sequence. The bound on x m given by the lemma follows from the Hadamard-Cauchy root 
formula for the radius of convergence of f(Z) around Z = 0. If K% is taken slightly greater than the 
inverse of the radius of convergence and K\ > 0, the bound \X m \ < KiK™ holds for large enough 
m. So K\ can be chosen to make the bound hold for every m = 0, 1, 2, . . . 

An explicit lower bound for the radius of convergence in terms of cq, . . . ,ci can be determined 
using the implicit function theorem proved by Lindelof using his majorant technique |13l p. 63] 
[20]. □ 

We are now ready to prove convergence of the formal psi series of Lemma 13.21 

Theorem 4.5. Consider the formal psi series of Lemma \3.2\ with rj = \og(b(t — to)) an d \b\ = 1- 
The branch cut is the segment 

{t - bp\p > 0}. 

Then the psi-series expansions for x(t), y(t), and z{t) given by (|1.2D or (13.11) converge uniformly 
and absolutely on the disc \t — to I < r with r > and with an open neighborhood of the branch 
cut excluded from the disc. In general, r will depend upon both C and D, which are the two 
undetermined constants in the psi series. 

Proof. We will give the proof for z(t). The proofs for x(t) and y{t) are similar. 

Excluding a neighborhood of the branch cut means that a neighborhood of to is excluded from 
the domain of convergence. Therefore R-2(t — to) 2 and R—i(t — to) are both bounded on the 
domain of convergence. The other reason for excluding a neighborhood of the branch cut is to 
ensure that rj is well-defined. 

By the definitions of \R m \ and \X m \ given at the beginning of this section, 

\Rm{r})(t - t ) m \ < \Rm\ max(l, |log6(t - t ) + C\ L(m+2)/2j ) |t - t | m 

< |X m |max(l,|log6(t-t ) + C| L(m+2)/2j ) |t-t | m 

< K x K?wbx.(i, |logb(t - t ) + c| L(m+2)/2j ) |t -to| m , 
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where m > and \b\ = 1. The first inequality above uses Lemma 13. II and the third inequality uses 
Lemma 14.41 

Choosing an r > such that 

r<l/K 2 and rfllog r\ + it + \C\) < 1/K 2 (4.8) 

is sufficient to ensure uniform and absolute convergence. The tt in (|4.8p is explained by the inequality 
[log 6(t — to) I < | log \ t — toll + 11 ■ A choice of r in accord with (|4.8p suffices for the convergence of 
the psi series for x(t) and y(t) as well. □ 

A further argument is required to show that the convergent psi series actually satisfy the dif- 
ferential equation, a point that seems to have been overlooked on occasion. When the psi series for 
x(t), y(t) and z(t) are substituted into the Lorenz system (jl.ip . the summation and multiplication 
of psi series on the right hand side is justified by standard results on rearrangements of absolutely 
convergent series. To justify the differentiation of psi series on the left hand side, we mention that 
the uniform convergence of a sequence of analytic functions on an open set implies the uniform 
convergence of the derivatives on any compact subset of that open set |28t Theorem 10.28]. We 
can now state the following theorem. 

Theorem 4.6. The psi series for x{t), y{t) and z(t) given by (jl.2p or (|3,ip . whose formal validity 
is asserted by Lemma \3.S\ satisfy the Lorenz system (jl.lh in the disc \t — 1$\ < r, with the branch 
cut excluded, for some r > 0. In general, r will depend upon both C and D, which are the two 
undetermined constants in the psi series. 

Levine and Tabor [19] raised the possibility that the locations of the singularities of an orbit of 
the Lorenz system may have accumulation points in the complex t-plane. Theorem 14.61 shows that 
psi series singularities cannot be accumulation points. 

So far, in results such as Lemma 13.21 and Theorem 14.61 we have regarded the psi series as 
functions of t. It is useful to consider them as functions of rj, where r\ = log6(t — to) gives a 
parametrization of the Riemann surface that gets rid of the branch cut in the t-plane. To be 
specific, we assume 9(to) < and b = —i. In that case, we have (t — to) = zexp(ry) and the psi 
series (|3.ip take on the form 

oo oo oo 

x (r,) = Y, i m P m (v)e m11 y(v) = im Qm{ri)e mv z(r,) = £ i m R m {r!)e m \ (4.9) 

m=— 1 m=— 2 m=—2 

with P m , Q m , R m being polynomials in which n always occurs in the group (rj + C). Every time t 
passes through the branch cut of log(— i(t — to)), n increases or decreases by 2-iri. Because rj and C 
always occur in the group (ry + C) in P m , Q m , R m , we can allow for other branches of log(— i(t — to)) 
in the psi series of (13. ip or (11.21) by keeping the principal branch of the logarithm in the definition 
of rj and incrementing C by an integer multiple of 2-Tri. The change in the estimate for the radius 
of convergence r of Theorem 14.61 for these other branches will then be in accord with (|4.8p (note 
that K 2 depends only on D). 

If the domain of convergence of the transformed psi series (|4.9p is considered in the 7/-plane, 
the choice of the principal branch of log(— i(t — to)) implies — ir < 9(r?) < tt and the r estimated 
by Theorem 14.61 implies 3ft(??) < logr. Thus the region of convergence of the principal branch will 
be a semi-infinite rectangle in the 77-plane. To pass to other branches, we keep C fixed and allow 



16 



Figure 4: Schematic plot of domain of convergence of the psi series in the rj plane, as implied by 
Theorem 14,51 The shape of the region is given approximately by 5J(r/) ^ — log|m| for integers m 
of large magnitude and — ir + 2irm < 9(r/) < ir + 2irm. 

the imaginary part of rj to be arbitrary. For r/ corresponding to different branches, one has to use 
different estimates for r as explained in the previous paragraph. Therefore the estimated domain 
of convergence of the transformed psi series (|4.9p will be a union of semi-infinite rectangles as in 
FigurelU If we start at the principal branch of log(— i(t— to)) and cross its branch cut m times, then 
by (|4.8jl r ~ 1/(^2 27r \m\) for large integers m. For such a branch — ir + 2irm < ^s(rj) < ir + 2irm 
and K(r/) ^ — log |m| for convergence, which gives an approximate idea of the shape of the domain 
sketched in Figure HI 

4.2 Remarks on theorems of Hille and Smith 

In [12] . Hille proved that the plane quadratic system 

dx/dt = x(ao + a\x + a2y) 
dy/dt = y(b + b\x + b 2 y) 

has a logarithmic psi series singularity if (a± — 6i)(o2 — 62) / («i^2 — 02^1) is a positive integer. Smith 
[30] generalized that result to plane polynomial systems. Smith's proof is based on a reduction to 
results proved early in the 20th century for Briot-Bouquet systems. These results are summarized 
in Sections 12.5 and 12.6 of Hille's book [13J. 

One difference between the results of Hille and Smith for plane polynomial systems and Theorem 
14.61 is as follows. The singular solutions for plane polynomial systems look like simple poles near 
the singular point. The singularities of the Lorenz system implied by Theorem 14.61 look like double 
poles. 

In [11], Hille proved the existence of logarithmic psi series solutions for the Emden- Fowler system 
d 2 y/dt 2 = i~ 2 /Py 1+2 /P for p = 2. At the end of the paper, Hille discussed the difficulty of extending 
his technique and noted remarks by a referee suggesting a proof of existence of logarithmic psi series 
solutions for positive integral 2p. Like Smith's proof for plane polynomial systems, the suggested 
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proof goes through a reduction to a Briot-Bouquet system, but no complete proofs are found in the 
literature as far as we are aware. The result for positive integral 2p was stated as Theorem 12.4.2 
in Hille's book |13| . Hille mentioned that "the various proofs are nasty," while omitting them. 

The proofs using reduction to Briot-Bouquet systems are difficult to follow in their entirety, 
partly because they depend so crucially on results proved long ago. It appears that use of the 
Laplace transform and the implicit function theorem will give simpler proofs for plane polynomial 
systems and complete proofs that are not so nasty in the case of the Emden-Fowler system with 
2p a positive integer. 

Theorem 4 of Smith's paper [30J states that all singularities of real solutions of certain plane 
polynomial systems must be of the form determined in Theorem 3 of that paper. The statement 
occurs again as Theorem 12.6.3 of |13j . Smith's proof begins with an ingenious change of variables. 
Near the end of the proof, we find the argument "in the case when A > , the arbitrary constant 
c in (20) can be chosen to fit this solution in the neighborhood of £ = 0." We are unable to 
follow that argument and believe it requires substantial explication at the very least. 

5 Complex singularities and the Lorenz attractor 

If t = to is a singularity of the Lorenz system (jl.ip . the solution must diverge to infinity as the 
singularity is approached. 

Theorem 5.1. Let 7 fe a Lipshitz curve in the complex t plane that approaches to at one of its 
two endpoints. Let (x(t),y(t), z{t)) be a solution of the Lorenz system (II. lj) defined for t € 7. If to 
is a singular point, then 



as t approaches to along the curve 7. 

Proof. Denote \x(t)\ + \y(t)\ + \z(t)\ by rp Consider the set of all complex (x,y,z) in the region 
\x — x(t) \ + \y — y{t) \ + \z — z(t)\ < b for some b > 0. Then the sum of the absolute values of the 
right hand sides of the Lorenz system (jl.ip is bounded by 



where 10 + 10 + 28 + 1 + 8/3 < 52 explains the first coefficient. 

Theorem 8.1, Chapter 1, of pQ (also see Theorem 2.3.1 of [13]) with a = 00 and M and b as 
above implies that the solution admits a unique analytic continuation to all t' in the disc \t' — t\ < R 
with 

p_ b 

52{r t + b) + 2(r t + b) 2 ' 

Taking b = r u we get R = 1/(104 + 8r t ). 

Being a singular point, to must lie outside the disc of analyticity. Therefore |to — t\ (104 + 8r t ) > 
1. Taking the limit t — > to along points on 7 completes the proof. 

The curve 7 is assumed to be Lipshitz to ensure uniqueness of the solution. □ 




(5.1) 



M = 52(r t + b) + 2(r t + b) 2 , 
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AB 



±i0.1714501006 



AAB 



±i0.1617621257 



AAAB 



±i0. 1563426260 



AABB 



±i 0.1636066901 



Table 2: Imaginary parts of the singular points closest to the real line in the i-plane. 

Theorem 15.11 proves that as the singular point to of the Lorenz system is approached, the 
magnitude of the solution must diverge at a rate that is at least as great as 0.125/ \t — to\. In fact, 
if the answer to Question 11.11 is yes and the singularities of the Lorenz system are all given by psi 
series of the form (|1.2|) . the divergence would be proportional to 1/ \t — to\ 2 . 

Theorem 15.11 is used to prove the theorem below. 

Theorem 5.2. Consider a trajectory of the Lorenz system (II, ip which is real for real values of t. 
In particular, assume that the state (x(0), y(0), z(0)) at t = is real. Then there is no singularity 
at any finite and real value oft and the solution is defined for all real values oft. 

Proof. Let Q = x 2 + y 2 + z 2 . Prom (jl.lj) , we have 



The matrix 1-norm of the symmetric form on the right hand side is bounded by 58 and so are the 
magnitudes of its eigenvalues. Therefore, \dQ/dt\ < 58Q and 



At this point it appears as if the proof can be completed using the Gronwall inequality (Theorem 
1.6.6 of [13]) to deduce that Q(t) < Q(0) exp(58 \t\). However, the bound on Q(t) holds only if we 
assume the existence of the solution, which is what we set out to prove. 

In circumstances such as these, oscillatory singularities for which the solution does not tend to 
a limit as the singular point is approached must be ruled out — an important point that goes back 
to Painleve |13l Chapter 3] . Theorem 15.11 forces the norm of the solution to diverge near a singular 
point thus making it possible to complete the proof. □ 

Theorem 15.21 is implied by Theorem 2.4, part (i), of [8]. In fact, Theorem 2.4 of [8] is a sharper 
result as it implies that Q(t) < C exp(20 \t\) for some constant C independent of t. We have given 
a proof that brings out the connection to the nature of the singular points. 

So far we know that the Lorenz system has singularities represented by logarithmic psi series 
and that the solution must diverge as a singularity is approached. But do solutions such as the one 
shown in Figure Q] have complex singularities and are they represented by psi series? 

Using numerical methods based on [26], we found the complex singularities closest to the real 
line of a few solutions listed in Table [2j Those solutions are all of course real for real t. They 
are assigned the labels AB, AAB, AAAB and AABB following the convention explained in the 
caption to Figure [TJ From Table [21 we see that the complex singularities are located at a distance 
greater than 0.037 from the real line, in agreement with Theorem 2.3 of Foais et al. [8]. In addition 
to computing the location of the singularities, we have verified numerically that their form matches 
the formal development of psi series given in Section 3. This numerical work will be described in 
detail elsewhere. 



dQ/dt = 2(-10j; 2 - y 2 - 8z 2 /3 + 38xy). 
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6 Conclusion 



Given that the Lorenz system (II. 1|) has resisted mathematical analysis on the real line, one may 
say that it is natural to think of t as a complex variable and analytic functions of t. When 

the solutions of the Lorenz system are viewed as analytic functions, it is natural to begin their 
investigation by looking at their singularities. We have given a complete formal development of 
singularities in the complex t-plane, proved convergence of the psi series representations using a 
new technique, and proved that the psi series indeed satisfy the Lorenz system. The development 
of the analytic theory appears to be a fascinating avenue for further investigations. 

Our suggestion that the mathematical analysis of the Lorenz system (jl.ip could be a problem 
in analytic function theory is an attempt to complete the circle, because the geometrical theory of 
differential equations, in which the Lorenz system is a famous example, sprang out of problems in 
analytic function theory — a fact that is not too well-known. More specifically, the stable manifold 
theorem, which is undoubtedly fundamental to the geometrical theory, was first proved to under- 
stand the solution of dz/dw = P(z, w)/Q(z, w) in a neighborhood of z = w = when P and Q are 
bivariate polynomials with P(0,0) = Q(0,0) = [13, p. 97] [271 1880]. 

The properties of analytic functions x(t) which satisfy the nonlinear Riccati equation dx/dt = 
fo(t) + fi{t) x + f2(t)x 2 , where the fi(t) are rational in t, is a well-studied topic. All the movable sin- 
gularities of the Riccati equation are poles and the dependence of its solution on the undetermined 
constant is given by a fractional linear transformation. For the Lorenz system some of the movable 
singularities have psi series representations of the form determined in Section 4. The dependence 
of these psi series solutions on the undetermined constants is much more complicated than for the 
Riccati equation. 

Another well-studied topic is the classification of second order nonlinear systems all of whose 
movable singularities are poles. The Painleve classification has been presented with lexicographic 
thoroughness by Ince [16] . There appear to be few classification results for third order systems such 
as the Lorenz system. Studying a specific system will probably sidestep many difficulties of the 
classification problem. In any event, the movable singularities of the Lorenz system are not poles. 
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